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Abstract 

We solve a unified integral equation to obtain the x, Qt and Q dependence of the 
gluon distribution of a proton in the small x regime; where x and Qt are the longitudinal 
momentum fraction and the transverse momentum of the gluon probed at a scale Q. 
The equation generates a gluon with a steep behaviour, with A ~ 0.5, and a Qt 
distribution which broadens as x decreases. We compare our solutions with, on the one 
hand, those that we obtain using the double-leading-logarithm approximation to Altarelli- 
Parisi evolution and, on the other hand, to those that we determine from the BFKL 
equation. 
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1. Introduction 

Understanding the details of the small x behaviour of parton distributions is one of the most 
challenging problems of perturbative QCD [|l], Moreover this topic has recently become of 
particular phenomenological interest with the advent of measurements of deep inelastic scat- 
tering at the high energy electron-proton collider, HERA, which have opened up the small x 
regime 0]. There now exist data for the proton structure function F2 for x as low as x ~ 10~^. 
These measurements reveal a significant rise of F2 with decreasing x, which has been taken to 
signal novel phenomena, although whether, in fact, this is the case or whether conventional 
explanations suffice, remains to be settled. 

First we recall the conventional treatment of "hard" scattering processes involving hadrons 
at moderate values of x, say x > 0.05. Then the observable quantities are calculated in pertur- 
bative QCD using the mass factorization theorem [Q in which the coUinear singularities, which 
occur in the partonic subprocesses, are absorbed into universal parton densities. To be specific 
let us take as an example the longitudinal structure function^] of the proton Fl{x, Q^) and, for 
simplicity, consider only the gluon^ partonic constituent. Then we have 



where we have chosen the mass factorization scale to be the "hard" scale of the process. The 
absorption of the collinear logarithmic singularities make the gluon density g{x, Q"^) "run" , with 
a dependence given by the Altarelli-Parisi evolution equations p . That is perturbative QCD 
does not determine the gluon absolutely but only its evolution from a non-perturbative input 
form. Altarelli-Parisi evolution resums the leading as\og{Q'^ / Ql) contributions. In a physical 
gauge the dglo^ {Q"^ / Ql) contribution can be associated with a space-hke chain of n gluon 
emissions in which the successive gluon transverse momenta are strongly ordered along the 
chain 0, that is k'^^ -C . . . -C k'^n ^ ■ The next-to-leading order contribution corresponds 
to the case when a pair of gluons are emitted without strong fc^-ordering (and iterations of this 
configuration). Then we have a power of as unaccompanied by logiQ"^ /Ql). A key feature of 
this conventional partonic approach is that the gluonic structure function, Fl (which at lowest 
order arises from the subprocess 7*(7 qq) is calculated assuming that the incoming gluon has 
negligible transverse momentum (and hence virtuality) as compared to the scale of the hard 
process. That is, on account of strong-ordering, we are able to work in terms of the density 
g{x, Q^) of the gluon integrated over its transverse momentum Qt- 

At sufficiently high electron-proton cm. energy, ^/s, we encounter a second large variable, 
\/x ~ s/Q'^^ and in this small x regime we must resum the leading as\og{l/ x) contributions. 

^We choose Fl as our example, rather than F2, because we are interested in hard scattering observables 
which are driven directly by the gluon. We could have used F2 but then we would have to deal with the 
collinear singularities of F2 associated with the g qq transition. These singularities are absorbed into the 
universal sea quark distributions. The simplification in which we neglect the sea quarks would therefore have 
been incomplete, even at leading order, and we would have had to broaden the discussion. 

^The gluon dominates the other partons in the small x regime, which is the main concern of our study. 
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The key ingredients of the QCD framework in this regime are the high-energy Qr-factorization 
theorem ^ and the BFKL equation ^ ^ for the gluon distribution, F{x,Qt), unintegrated 
over its transverse momentum Qt- Using Qr-factorization our sample observable is now given 

by 

^1 dx' f (fQr , ^ fx Qj 

where Fl denotes the gluonic structure function calculated using an off-shell {Q"^ 7^ 0) gluon. 
If we were to return to the strongly-ordered transverse momentum configuration then 



FL-^FL = FL{x/x',0,as{Q^)) (3) 

and 

xg{x, Q') = [ ^F{x, QtMQ - Qt), (4) 

though, of course, it is inappropriate to work in terms of the familiar integrated distribution 
g{x, Q"^) at small x. 

The BFKL equation for the unintegrated gluon density F{x, Qt) sums the leading as log(l/x) 
contributions. The strong-ordering in transverse momenta is no longer applicable and we now 
have a "random walk" or diffusion in /c^ as we proceed along the chain. The enlarged fcy phase 
space leads to a x~^ growth, with decreasing x, where A ~ 0.5. Indeed the observed behaviour 
of F2 appears consistent with the precocious onset of this leading log(l/a;) behaviour of the 
gluon [0, but the definitive confirmation that this is the case must await the calculation of 



the subleading corrections to the BFKL equation. Clearly the BFKL equation, which resums 
the leading log(l/x) contributions, has a limited region of validity. In principle, it is restricted 
to the region a5log(l/x) ~ 0{1) and 0:5 log(Q^/(5o) ^ where Qq indicates the boundary 
of the non-perturbative domain, Qg ~ 1 GeV^. To make progress we need to know how the 
BFKL formalism links with the conventional Altarelli-Parisi dynamics at larger x and large 

log(QVQg)- 

A theoretical framework which gives a unified treatment throughout the x, Q^ kinematic 
region has been provided by Catani, Ciafaloni, Fiorani and Marchesini [11-15]. The resulting 
equation, which we shall call the CCFM equation, treats both the small and large x regions 
in a unified way. The equation is based on the coherent radiation of gluons, which leads to an 
angular ordering of the gluon emissions along the chain. In the leading log(l/x) approximation 
the CCFM equation reduces to the BFKL equation, whereas at moderate x the angular order- 
ing becomes an ordering in the gluon transverse momenta and the CCFM equation becomes 
equivalent to standard Altarelli-Parisi evolution |Q. The angular ordering introduces an addi- 
tional scale (which turns out to be essentially the hard scale Q of the probe), which is needed 
to specify the maximum angle of gluon emission. Thus we must work with a scale-dependent, 
unintegrated gluon density F{x, Qt, Q)- At very small x the angular ordering does not provide 
any constraint on the transverse momenta along the chain and F becomes the Q-independent 
gluon of the BFKL equation. 
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The aim of our paper is to study the CCFM equation in detail and, in particular, to 
obtain numerical solutions so as to reveal the small x behaviour of the gluon distribution in 
the proton. In Section 2 we sketch the derivation of the CCFM equation, while in Section 
3 we study its properties in the important small x region. We make approximations which 
both simplify the discussion and facilitate the solution of the equation. To gain insight we 
study both the "folded" and "unfolded" versions of the equation. In Section 4 we present 
numerical solutions F{x, Qt, Q) of the CCFM equation and compare them with solutions that 
we obtain from solving the two limiting versions of the equation. That is we compare the 
CCFM solutions with those obtained from (i) the double-leading-logarithm approximation to 
the Altarelli-Parisi equation and (ii) the BFKL equation. Section 5 contains our conclusions 
on the small x behaviour of the gluon. 

2. Coherent branching : the master equation for the gluon 

The perturbative evaluation of physical QCD quantities in general, and parton distributions 
in particular, is complicated by the presence of large logarithms which arise from the emission 
of both soft and collinear gluons. The origin of the large logarithms can be seen from Fig. 1. 
The differential probability for emitting a gluon of 4-momentum q is of the form 

dP ~ a A# (5) 

Zg Qrp 

where qt is the transverse momentum and Zg is the longitudinal momentum expressed as a 
fraction of the momentum of the parent gluon. 

Here we focus attention on the gluon distribution within a proton. To predict the correct 
behaviour of the distribution it is necessary to resum the large logarithms which arise not just 
from single but from multigluon emissions to all orders in a^. A typical contribution is shown 
in Fig. 2, where a gluon of low space-like virtuality evolves to higher virtuality and lower energy 
by successive gluon emission. It can be shown that the emissions are coherent in the sense that 
there is angular ordering, 6i > 6i-i, along the chain, where di is the angle that the i*'* gluon 
makes to the original direction [11-16]. Outside this region there is destructive interference such 
that the multigluon contributions vanish to leading order. We speak of coherent branching. 

As mentioned in Section 1, due to the presence of angular (rather than transverse momen- 
tum) ordering of the emitted gluons we need to expose the transverse momentum, Qt, depen- 
dence of the probed gluon. That is we work in terms of the scale {Q) dependent "unintegrated" 
gluon density F{x, Qt, Q), which specifies the chance of finding a gluon with longitudinal mo- 
mentum fraction x and transverse momentum of magnitude Qt [|13- The integral equation for 
F{x, Qt, Q), which effects the summation of the large logarithms, can be approximated on the 
one hand, to yield the BFKL equation at small x, where F becomes independent of Q and, on 
the other hand, to yield at moderate x the Altarelli-Parisi (or CLAP) evolution equation for 
the integrated distribution g{x,Q'^). It is this master equation which we wish to investigate 
and to solve. 
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It is necessary to outline the derivation of the equation |11, 12, 13 1. A crucial development 
has been the proof of the soft gluon factorization theorems which allow the inclusion, not 
only real gluon emission, but also the virtual emission contributions which tame the singular 
behaviour in the two boundary regions z —>■ and 2; — > 1. A recurrence relation can then be 
obtained which expresses the contribution from n-gluon emission in terms of that from n — 1 
emission. In essence, we integrate the differential probability dP for the emission of the extra 
gluon over the relevant region of phase space, where (H) takes the explicit form 

dP = AsPdz^Q{e-e'). (6) 

It 

Here — 6') reflects the angular ordering and P is the gluon-gluon splitting function 



P = as 



^ +AMs--'2 + zil-z] 



1 — z z 



(7) 



We define 0.3 = Cac^sIt^ = 3as/7r. The multiplicative factors and Aats, known as the 
Sudakov and non-Sudakov form factors, arise from the resummation of the virtual corrections. 
They cancel the singularities manifest as 2; — 1 and 2; — respectively. These form factors 
have exponential forms which we present in eqs. (|10D and ([TT|) below. 

We use to obtain a recursion relation expressing the distribution JF„ in terms of Tn-\- 
The distribution Qt-, z, q) corresponds to the n-rung ladder diagram, where the variables 

are defined in Fig. 3. Following refs. P3| , [T^ we impose the angular-ordered constraint by 
introducing rescaled transverse momenta 

q = ^ 9E', q' = ^ e'E" (8) 

^ l-z ' ^ l-z' ^ ' 

where 1 — 2; is the longitudinal momentum fraction of the gluon emitted at angle 9 and E' is the 
energy component of the exchanged gluon with spacelike momentum x'p. Here we have used 
the small angle approximation, tan^ 0. The coherence constraint 9 > 9' therefore implies 
q > z'q' and so (|D gives 



Tn{x, Qt. z, q) = f dz' f ^e(g - z'q')As{q, z'q')P{z, q, Qt)^„-i (-, Q't, z' , q') , (9) 
Jx/z J Tiq'^ \z J 



see Fig. 3. Due to the presence of angular ordering, it appears that we have to consider a less 
inclusive structure function than the unintegrated distribution Qt, Q) itself. To be precise 
we have exposed not only the x, Qt of the probed gluon, but also the z, q dependence which 
specifies the previous emitted gluon. Note that the variable in Tn-i is the magnitude of the 
vector sum + {1 ~ z)q and, so, in principle, the angular integration in d'^q' is non-trivial. 

The Sudakov form factor is given by 

A^(g, z'q') = exp f- ^ /' dx-^] . (10) 

\ J(z'q')^ Jo 1 — Z j 
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The region of integration corresponds to the angular-ordered region from the angle 9' of emission 
of the (n — 1)*'^ gluon to 9 of the n^^ gluon. Indeed A5 can be interpreted as the probability for 
not emitting a gluon in this angular region. This observation is consistent with which gives 
the differential probability for single gluon emission in an element of phase space specified by 
q^] the factor A5 ensuring that there is no prior emission. 

The splitting function P{z,q,QT) is given by (|^) in which the 1/z singularity is screened 
by virtual corrections contained in the non-Sudakov form factor 

CrZQ rl 2' r dh'^ \ 
-asj^ —J^QiQl-e)Qik-z'q)j (11) 



where 

r 1 if {Qr/q) > 1 

zo = I Qr/q a z < {Qr/q) < 1 

[ z if {Qr/q) < z. 

Unlike A5, the non-Sudakov form factor Ajvs is not just a function of the branching variables, 
but depends on the history of the cascade via 

QT = \qT + q'T + qT + ■■■V (13) 

Actually the recursion relation (j^) is satisfied by a more inclusive distribution Qt, Q) 
in which the 2;, q dependence of JF„ is integrated over, subject to the maximum angle specified 
by Q |T^. That is, if we introduce a distribution for n-gluon emission defined by 

Fnix, Qt, Q)= f dz ( -^QiQ - zq)AsiQ, zq)J^nix, Qt, z, g), (14) 

then the recursion relation becomes 

72 



Fn{x,QT,Q)= / — -^e{Q-zq)As{Q,zq)P{z,q,QT)Fn-ii-,Q'T,q). (15) 
Jx z J Tiq \z J 



Ix z J nq'^ 

Finally we obtain the (scale-dependent unintegrated) gluon density by summing over all gluon 
emissions 

00 

F{x, Qt, Q) = E F'nix, Qt, Q). (16) 

n=0 

From ([l5[) we find 

F{x,Qt,Q) = F\x,Qt,Q) + 

+ [dz J ^Q{Q - zq)As{Q,zq)P{z,q,QT)F(^^,Q'T,q'j, (17) 
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with Q'rp = IQj, + (1 — z)q\. The inhomogeneous or "no-rung" contribution, F^, may be 
regarded as the non-perturbative driving term. This basic integral equation (p!7[), for the gluon 
structure function F{x, Qt, Q), which we have called the CCFM equation after its originators, 
is the starting point of our analysis. It may be approximated both at moderate x to yield the 
Altarelli-Parisi evolution equation, and at small x to yield the BFKL equation, as we shall now 
show. 



Qt, Q) = F\x, Qt, Q)+ols [' — f ^Q{Q-zq)^Ns{z, g, Qt)F f-, |Qr + (18) 

J X Z J TXQ \Z / 



3. Approximations in the small x region 

To explore the structure of the gluon in the small x region we may approximate (|T^) by 

_ dz f d'^q ^ 
Txq^ 

where we have set A 5 = 1 and retained only the 1/ z term in the splitting function P. We have 
also approximated (1 — z)g by g in the argument Q'rp of F. In this small x limit, the variable q 
reduces to the transverse momentum of the emitted gluon, see eq. (H). This is the CCFM 
integral equation [|ri|, |l^] which we solve to find the x, Qt and Q dependence of the gluon 
distribution. The procedure that we adopt and the results that we obtain are presented in 
Section 4. 

Also in Section 4 we compare the results obtained from the CCFM equation with those 



obtained in the double leading logarithm (DLL) approximation in which (p!8| ) reduces to 

F{x, Qt, Q) = F%x, Qt, Q) + ast - I ^e(Q - q)F + q\,q). (19) 

Jx z J irq"^ \z J 

To be precise the DLL approximation is obtained by setting A^vs = 1 and by replacing the 
angular ordering, Q{Q — zq)^ by ordering in transverse momentum, Q{<Q — q). This procedure 
becomes equivalent to the conventional DLL approximation for xg{x,Q'^^ after we integrate 
over Qt, as in (^). Here we have the advantage that we can also display the Qt dependence of 
the gluon distribution. 

Before we present our numerical results it is informative to gain insight into the structure 
of the CCFM equation, ([TB|). We can simplify the discussion by rewriting ( ]TB| ) in terms of the 
moment function 

F^{(Qt, Q) = t x'^-^Fix, Qt, Q). (20) 

JO 



We obtain 



where 



dj^ 

FUQT,Q) = F!i{QT,Q) + as J :^H^{Q,QT,q)FU\QT + qlq) (21) 



H^{Q,QT,q) = C dzz^-^Q{Q-zq)/\Ns{z,q,QT) (22) 

JO 

= Q{Q-q) C dzz-'^ANs{z,q,QT) 
Jo 

rQ/i 

+ e{q-Q) dzz^-^^Ns{z,q,QT). (23) 
Jo 



The curtailment of the range of integration in the second term is seen to be a direct consequence 
of angular ordering. 



(a) Unfolding the equation 

It is useful to unfold the kernel of ( pT]) so that the real emission and virtual corrections 
terms appear on equal footing, i.e. to the same order in as- On the one hand this will allow 
us to make correspondence with the BFKL equation, and on the other hand it will guide us 
to the correct structure of the driving term F^{x,Qt,Q) to be used as input for the CCFM 
To unfold the kernel we first integrate 

1 



equation, (18 



UJ 



e{Q-q) + {Q/qr&{q-Q)A^s{z 



3D by parts 

Q/q,q,QT 

-— fdz I %e{Q - zq)e{k - zq)Q{Ql - e)Ar,s{z, g, Qt) 
UJ Jo J 

- \<d{Q -q) + {Q/qT - Q)Ans {z = Q/q, q, Qt 

- ^ / ^0(Qt - k')H^{mm{Q, k}, q, Qt) 

31). We insert (M) in 



where the last term has been simplified using 
term in terms of F^, — F^,. We obtain 



(24) 

If) and express the final 



F.(Qt, Q) = F^{Qt, Q) + — f %Q{Ql - q')F^iQT, min{Q, q}) 

UJ J q"^ 



+ — 



as f d'^q 



UJ J rrq^ 



<d{Q-q) + {Q/qYQiq - Q)Ans{z = Q/q, g, Qt) 



- e(Q|-g2)F^(Q^,min{Q,g}) 



F^{\QT + qlq) 
(25) 



It remains to unfold the non-Sudakov form factor Ajvs- In the Appendix A, we show that (|25D 
then becomes 



UJ UJ J nq"^ 

as f d^q 



Q{Q-q)+(^] e{q-Q) 



FUlQx + Qlq) 



-- mml-q')F.{QT,q) 
UJ J Trg^ 



+ —Q{Qt-Q / ^ —q 
UJ JQ^ Tiq^ 



Q'rd'q ,dFUQT,q), fq' 



A 2 I ^2 

Q2 Tiq^ dq^ \Qt, 



'q: 
.1. 



(26) 



where the driving term of the unfolded equation is related to that of the folded equation by 

1 ^n.„ „n.„ as fd^q^ 

UJ " ' ' ■ UJ J vrg^ 



-Fl{iQT,Q) = F^{Qt,Q) + ^ I ^e{Q'T-q')F^{QT,q) 



«5 
UJ 



Ql d'q^,dF^iQT,q)_^^^f q 



Q2 vrg" 
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dq"^ 



Ql 



'q; 
.1. 



(27) 



The second term on the right hand side of ( p6D corresponds to real gluon emission without 
any "non-Sudakov" damping, but with angular ordering included (cf. (|23|)). The following two 
terms are the unfolded virtual corrections, which, if resummed, would lead to the non-Sudakov 
form factor A^s- The inhomogeneous term in the unfolded equation, (PU]), should not 
contain any virtual corrections, but in the folded equation, ([TsD , may. 

To solve (|27|) for F^ in terms of F^ it is simplest to use the following representation for the 
driving term of the unfolded equation (|26|) 

-F^{Qt, Q) = f dx f'^ f^e{Q- (-, \q + Q^\,q) . (28) 

u! Jo Jx z J T:q^ \z J 

Then the solution of (0) is 

^°(Qt, Q) = t dx x-~' C- I ^e(Q - qz)Ar,s{z, q, Qt)$° (-, \q + QtU) • (29) 
Jo Jx z J nq^ \z J 

The above representation, (|28|), mirrors relation (14) which, for n = 0, expresses the driving 
term F^ in terms of the non-perturbative input Tn=o- In section 4 we will use (^) to specify 
the driving term F°(a;, Qr, Q) of the (folded) CCFM equation ([T8|) in terms of an assumed form 
for the non-perturbative input function 

(b) The BFKL limit 

We may take the leading logl/x approximation of the unfolded equation, (p^), for the 
moments of the gluon distribution. This corresponds to retaining the leading terms in the 
a; — >■ limit. That is we set {Q/q)^ = 1 in (|26|), which then reduces to 

F.{Qt, Q) = -F°^(gT, Q) + ^j—, [f.{\Qt + qU) - F^{Qt, q)Q{Ql - (30) 
UJ UJ J Tiq '- -' 

where F^^ is the BFKL limit of F^ of (0). We see that the equation ( pUD generates a moment 
function F^^ which is independent of Q. If we transform back from moment (uj) space to x 
space, then we obtain the BFKL equation 

(91og(l/x) Tcq'^ 

where we have neglected the derivative of the inhomogeneous term with respect to log(l/x). 
The two terms in the integral on the right hand side correspond to real and virtual gluon 
emission respectively. 

(c) Cancellation of the real and virtual singularities 

We first expose the cancellation of the singularities in the BFKL limit 0. To do this we 
rewrite ( pUD with the "unresolved" real emissions (i.e. emissions with < yU^) separated out 

as f d'^q 



J [f{x, \Qt + q\)- Q{Ql - q')F{x, Qt)] , (31) 



fuQt) = -f'^\Qt) + — I—, [f.{\Qt + qm^^' - q') - f^{Qt)Q{qI - q') 

UJ UJ J Tiq'^ L 

+ — l^,U\QT + qmq'-^% (32) 
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uo J T:q^ 



The singularities as /i^ ^ occur in the second term on the right hand side. To exphcitly show 
the cancellation between the real and virtual contributions, we approximate the "unresolved" 
real emission contribution using F^{\Qj, + qr|) ~ F^{Qt), which is valid at small q, and obtain 



as 



Ju? 



-F^{Qt)uJ{QI,^?) (33) 



where 



TD{Ql,^^'') = as\og{Ql/^^^). (34) 

The result (|33D is the residual virtual contribution to ujF^^{Qt) which remains after the cancel- 
lation of the real and virtual singularities. 



If we substitute (133|) into (|32D , we obtain 

uo + ud u + uJ J Tcq 
which, when we transform back to x space, becomes 



F.iQr) = ^^^^ + ^ / - ^^')F^i\QT + (35) 



F{x,Qt)= f-z^F'^(^-,Q^ + asf-z''^^Q{q'-^^')F(^-^Q^ + q\). (36) 
Jx z \z J Jx z J vrg^ \z J 

Equation (^) clearly corresponds to the "folded" BFKL equation in which we have resummed 
all the "unresolved" real emissions and all the virtual corrections. 

We recognise the Regge-like form z'^ of the non-Sudakov form factor A^^ in the folded 
BFKL equation, which screens the 2; ^ singularities. We can also make the identification of 
the form factor using (|llD, which in this case becomes 

= expj -a5log(l/z)log(Q|//x^) 

= z"^, (37) 

where uJ is given in (0). Note that the driving term F^^/uj of the unfolded BFKL equation 
( |32| ) is free of /i^ ^ singularities, whereas the driving term of the folded version of equation 
(|35|) or (|36|) contains, through ZJ, these //^ singularities. 

In the same way we can introduce a resolution /z^ into the unfolded CCFM equation, (pBD , 
so that 
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F^{Qt,Q) 



UJ 



as f d^Q 



e{Q-q) + i^ e{q-Q) 



eiq'-ix')FU\QT + q\,q) 



- mmi-emq'-p^')F.{QT,q) 

UJ J Tiq^ 



+ — Qt - <5 / , -^q — w-0 log -Tw 



? - 



(38) 



After resumming all the virtual corrections and "unresolved" radiation we obtain, in analogy 
with (M 



;2 

Tiq 



2' 



F^{Qt, Q) = FliiQr, Q; /i') + as I ^H^iQ, Qt, q; fi')F^{\QT + q\,q)Q{q' - 



with 



Hu;{Q,QT,q]fJ''^ 



dz z^-'Q{Q - qz)ANs{z, g, Qt; I^') 



(39) 
(40) 



where now the non-Sudakov form factor is of the form 



/ /■! rl7' r ftk"^ \ 
A;v5(^, g, Qt; /i') = exp {-as j^—j ^©(Ql - ^')©(fc - z'q)Q{k^ - ^?)\ . (41) 



The driving term of (^) has a similar form to ( ppj ) except that the introduction of the resolution 
cut-off leads to a factor 6(g^ — /i^) and to /S.ms being replaced by Aats of ( PI) 



^°(QT,Q;/i') = r^xx^-i / ^e(Q-gz)e(g2V)A;v5(^,g,QT;/^')$°(-,|q+QrU)- 

JO Jx z J Txq^^ \z J 

(42) 

In fact we do not need to introduce a lower cut-off [i^ on the dq^ integration in (p!8D , since 
the equation remains finite as g — 0. However, here we have demonstrated the modifications 
necessary if, for pragmatic reasons, a cut-off is introduced. In particular we see Ans of ([TT|) 
must be replaced by ^ns of (^Tj). The results will, of course, be independent of the choice of 
the resolution /i^, up to 0{^'^ /Q\). 

Although the unfolded CCFM equation, (^Bf) , reduces to the BFKL equation in the leading 
log(l/x) (or leading as/uo) approximation and although both equations are free from singular- 
ities at q^ = 0, the details of the cancellations are different |jl5|. In the BFKL case the real 
emission terms and the virtual corrections are individually divergent and have to be separately 
regulated by the q^ = //^ cut-off. On the other hand, for the CCFM equation both terms are 
finite, yet they generate additional powers of l/u. To be precise, when we solve eq. (|26| ) 
iteratively we find after n iterations that 



Fi''\QT,q) 



~— g"log"(g/QT) 



(43) 
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as g — > 0. The factor q'^ regulates the integrals at = and keeps both the real and virtual 
terms finite, yet they separately contain double logarithmic factors, (log^(l/x))", arising from 
the behaviour of the integral 

- /^g-log"f^U^ (44) 

The (log(l/a;))^" behaviour exactly cancels between the real emission and the virtual correction 
terms. 



4. Numerical solution of the CCFM equation 

We explore the structure of the gluon distribution F{x, Qt, Q) at small x by numerically 
solving the CCFM equation (0). To be precise we solve the equation in the presence of the fj,"^ 
resolution cut-off (cf. (p9D) 



F{x, Qt, Q) = F\x, Qt, Q; /i')+ / - / -4^se{Q-zq)e{q'-fi')A^s{z, q, Qt, t^')F{-, Q't, q) 

Jx z J q'^ z 

(45) 

with Q'j. = \Qt + (1 — z)q\, and where the non-Sudakov form factor A^s is evaluated using 
(Pip. We take the argument of as to be Qt since this value is usually assumed for small x 
studies involving the BFKL equation. 

From (^) we see that the driving term is given by 



F%x,QT,Q;fi')= C- I ^Q{Q-qz)Q{q'-fi')A^s{z,q,QT;fi')<^'(-,\q + QT\,q). (46) 
Jx z J Tiq^ \z J 



In general $° will contain some smearing in q. However, this is an inessential complication and 
it is sufficient to assume strong-ordering in transverse momenta and to take the q dependence of 
$° to be 5(q-|-Qy)- The xj z and Qt dependence of $° are chosen so that if the functions and 
Aj-^s were to be set equal to 1 in (^6]) , then F° would reduce to 3(1 — a:)^A^exp(— Q|,/Qq). The 
normalisation iV is fixed so that the gluon, integrated over all Q\, carries half the momentum 
of the proton. Q\ is taken to be 1 GeV^. These assumptions are equivalent to the choice 

/■I Ay 

F\x,Qt,Q;[i^) = Nexp{-Ql/Ql) _e(Q - QT^)e(Q| - /x^) 

Jx Z 

A (7 n n ■ ;;2N ^[3(l-a^A)'] /.yX 

alog{z/x) 

We solve ( ^5] ) by iteration starting from the input form given in (^7]). We restrict the 
iterative procedure to the domain^ Q'^,Q'^ > Q^ = 1 GeV^. We take the upper limit cut-off, 

^Strictly speaking the non-Sudakov form factor, (^l]), is also dependent on Qq to ensure complete cancellation 
of the real and virtual emissions. 
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Qf, on the integrations to be in the region 10^ - 10^ GeV^ though the results are insensitive 
to variations around and above these values. We use a lower cut-off /i^, although we saw in 
Section 3 that the CCFM equation is well behaved in the /i^ — limit. 

Clearly, therefore, the results should be independent of the choice of the resolution /x^, up to 
contributions of Oi^ji^ /Q'^) and Oi^ji^ /Q"^). This is well demonstrated by the sample of results 
shown in Figs. 4 and 5 which correspond to the choices /i^ = 10~^, 10~^, 0.5 and, 1 GeV^. The 
convergence to a stable result with decreasing ii^ is a test of the consistency of the solution. 
From now on we show results for the choice jj?' = 10"^ GeV^. 

(a) CCFM solutions compared with those of the DLL approximation 

Starting from the input of (^Tj) we iterate ( ^Sj) until a stable result for F{x, Qt, Q) is 
obtained. The numerical procedure is outlined in appendix B. We take Q\ = 10^ GeV^. The 
whole calculation is repeated using the double-leading-logarithm (DLL) approximation in which 
we replace angular ordering, Q{Q — zq), by ordering in transverse momenta, Q{Q — q), and in 
which we set A^vs = 1- 

Figs. 6, 7 and 8 respectively show a representative sample of results for the Q, Qt and x 
dependence of the gluon distribution in the small x regime. The three plots of Fig. 6 compare 
the dependence of F{x, Qt, Q) obtained from the CCFM and DLL equations, for Q"^ = 1, 10 
and 100 GeV^ respectively. We see that F obtained from the CCFM equation is less dependent 
on Q^ than the DLL values of F. In the DLL case the suppression of F in the region Q^ < Q^ 
is simply a reflection of the ordering of transverse momenta that is embodied in the equation; 
an ordering which is not, in fact, appropriate in the small x domain. Recall that in the BFKL 
leading log(l/x) limit F becomes independent of Q^. The CCFM solutions in Fig. 6 exhibit this 
behaviour for the larger values of Q^, but for smaller Q^ they show that non-leading log(l/a:) 
contributions begin to become important. The low Q^ region is where the physically motivated 
angular ordering embodied in the CCFM equation (but not in the BFKL equation) provides 
more of a constraint. 

Fig. 7 shows the Qt distributions of the gluon obtained from solving the CCFM equation and 
from the approximate DLL equation. As before the DLL solutions satisfy Q'^ ^ Q^ even at the 
smallest values of x, which again reflect the transverse momentum ordering, Q{Q — q), contained 
in equation (|19|). On the other hand the CCFM solutions become significantly broader in Qt, 



with decreasing x, on account of the more appropriate angular ordering constraint Q{Q — zq). 
The extensive Qt tail is a key property which characterises the gluon distribution in the small 
X domain. 

Fig. 8 shows the behaviour of the integrated gluon distribution, xg{x, Q"^) of (^) with lower 
limit Qq, as a function of x for various values of Q^. Note that the gluon distributions are 
generated radiatively from an input which is "fiat" at small x, (^61) , and so the rapid rise of 
xg with decreasing x (shown as continuous curves) is generated by the CCFM equation. To 
quantify the increase, we show in Fig. 9 the effective value of A, defined by 

xg{x,Q^) = Ax-^. (48) 
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For small x we see that the solutions converge to a typical x~^'^ behaviour, approximately 
independent of Q^, which, as we shall see below, is consistent with that obtained from the 
solution of the (leading log(l/x)) BFKL equation. 

The dashed curves in Figs. 8 and 9 show the characteristic double-leading-logarithm (DLL) 
small X behaviour 



xg{x, Q ) ~ exp 



2{e(Q^Q^)log(l/a;)}^ 



(49) 



appropriate to a "flat" input, that is xg{x, Ql) — > constant as x — * 0. In (^) we have omitted 
slowly varying functions of the argument of the exponential. We see that the steepness or 
"effective slope" A increases with via the "evolution length" 



aQ',Ql)= r^-^asiq'). (50) 



This behaviour is clearly evident in the DLL results of Figs. 8 and 9. The difference with 
the CCFM predictions for small x, a; ^ 10~^, shows the importance of implementing angular 
ordering in this domain. 

(b) Comparison with solutions of the BFKL equation 

To investigate the small x limit we compare the CCFM solutions with those that we obtained 
by solving the BFKL equation with the same driving term. To be precise we solve the unfolded 
BFKL equation 

F{x,Qt) = F'''^{x,Qt) + 

+ -^kwi — m^i — ^(iQ^TQpr''^ 



z 



with the driving term 

F°^(x,Qt) = 3(l-x)^Arexp(-Q|/Q2). (52) 

The integration region is restricted to > Qo = ^ GeV^. For completeness, we show in 
Appendix C that, for fixed coupling 0:5 and Ql = 0, the BFKL equation (^) can be rewritten 
in the form 



F{x,QT)=F'\x,QT)+as t- 1^ 

Jx z J Tcq'^ 



2, 

2 „2n 



F{x,\QT + q\)-e{Q'T-q')Fix,QT) 



(53) 



which we obtained previously, see (|3lD . 

Fig. 10 compares the properties of the integrated gluons obtained from solving the CCFM 
and BFKL equations, while Fig. 9 shows the resulting effective slopes A. Since the solution, 
F{x,Qt), of the BFKL equation is independent of Q, the Q'^ dependence observed for xg 
comes entirely from the Qt integration of (^. On the other hand the solution F{x, Qt, Q) of 
the CCFM equation has an intrinsic Q dependence arising from angular-ordering, Q{Q — qz). 
As a consequence we see from Fig. 10 that the CCFM gluon evolves faster with Q^. From Fig. 
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9 we note that the effective slopes A of the integrated CCFM and BFKL gluons are remarkably 
similar at small x. We conclude that the next-to-leading log(l/a;) effects included in the CCFM 
formalism have a comparatively weak effect on the x"^ behaviour, although we note that the 
onset of the form is more delayed for the CCFM solution. 

5. Conclusions 

We have solved a unified equation for the unintegrated gluon distribution which incorporates 
BFKL dynamics at small x and Altarelli-Parisi evolution at larger x. We called it the CCFM 
equation after its originators — Catani, Ciafaloni, Fiorani and Marchesini. Starting from a 
driving term based on a "flat" 3(1 - xf gluon with a narrow Qt distribution, exp( — Q't/Q'q)-! 
we used an iterative procedure to find the x, Qt and Q dependence of the gluon. 

We concentrated on the behaviour of the gluon in the small x regime. The key ingredients 
of the CCFM equation are the angular-ordering of gluon emissions and the presence of a non- 
Sudakov form factor. We found that the CCFM equation generates a gluon Qt, Q) with 
a singular x~^ behaviour, with A ~ 0.5, and a Qt distribution which broadens and develops a 
significant tail as x decreases. Moreover the angular-ordering introduces a dependence of the 
unintegrated gluon on the scale Q, especially at the lower values of Q^. Sample results are 
shown in Figs. 6-9. It is convenient to display the x dependence of the integrated distribution, 
xg of (H), although we should recall that the physically relevant quantity at small x is the 
unintegrated distribution Qt, Q)- 

We compared the CCFM solutions with the conventional DLL approximation in which 
angular-ordering is replaced by strong-ordering in the gluon transverse momenta and in which 
the non-Sudakov form factor is omitted, ^ns = 1- The gluon is then found to be much less 
steep with decreasing x and to have a narrower Qt distribution. We found that the DLL 
approximation starts to differ from the CCFM results in the region x ^ 10~^. 

We then compared the CCFM solutions with the solutions of the BFKL approximation, 
based on an equation in which the angular-ordering is ignored, and which therefore has unin- 
tegrated solutions F{x,Qt) which do not depend on the scale Q. In fact the gluon g{x,Q'^) 
obtained from the BFKL equation acquires its Q^ dependence entirely from the integration in 
(^). As a consequence we find that the integrated BFKL solutions evolve more slowly in Q^ 
than those obtained from the CCFM equation, which have an additional intrinsic Q^ depen- 
dence as shown, for example, in Fig. 6. Fig. 9 quantifies the x""*" agreement between the unified 
CCFM solution and the approximate BFKL solution. The agreement is remarkably good at 
small X, especially at the larger values of Q^. Both the CCFM and BFKL solutions have a 
behaviour xg ~ at small x, where the value of A is in the region of 0.5 with only a modest 
dependence on Q^, in contrast to the dependence of A on the evolution length for the DLL 
approximation, see Fig. 9. 

It is appropriate to add a word of caution. Our study has concentrated on obtaining the 
perturbative QCD predictions of the behaviour of the gluon at small x. In particular we have 
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made the small x approximation for the splitting function (|^) 

P ~ asA^s- (54) 

and we have set the Sudakov form factor A5 = 1. 

In other words we have neglected the singularity at z = 1 and its Sudakov suppression. 
Moreover we have neglected the effect of the flavour singlet quarks on the evolution of the gluon. 
Although these may be viewed as moderate-to-large x effects, they may, to some extent, feed 
through to small x leading to possible next-to-leading log(l/a;) contributions to F{x,Qt,Q)- 
Therefore, at this stage, our results should be regarded as illustrative of the main properties of 
the gluon distribution at small x, and as an indication of the respective domains of validity of 
BFKL and conventional Altarelli-Parisi dynamics. To obtain a quantitative prediction of the 
gluon for all x we must include a proper treatment of the z = 1 behaviour and include the 
quark distributions . 
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Appendix A : Unfolding the CCFM equation 

Here we derive the "unfolded" version of the CCFM equation, which is shown in (|26|) . We 
start by unfolding the non-Sudakov form factor A7V5 in (^51) . The relevant term is 

r ^ (-] Ar,s{z = Q/q, q, Qt)F^{\Qt + qU) 

'-'^ d?q (Q^ 



q''guq,qt) + r 



Ttq- \ q 



F^{\QT + qlq) 



with 



G.(Q, Qt) ^ (1) [A^^(^ = Q/q, g, Q^) - 1]FU\Qt + ^U), 



(55) 



(56) 



?2 vrg' vg" 

where we have simply subtracted and added 1 to A^s- For Q^ < Q^ the unfolding is trivial 
since, then Aj^si^: = Q / q^q^Qr) = 1, see (p!2D. This leaves the case Q'^ > . From the form 
of Ajsss we see that G^^Qt-, Qt) = 0, and so we may write 



G.(Q,Qt) 



12 



Ql dQ 
15 



OG^iQ'^Qr) 



d log Q 



12 



(57) 



If we differentiate 



we obtain 

- (Pq 1 d/\Ns{z = Q'/q,q,QT) 



d log Q 



12 



2 r,i^ 



Q'2 vrg" q 



d\ogQ' 



FU\QT + q\,q) 



y2 



1 1^" 



(58) 



91ogQ'2 (91ogQ'2 

where the last equality is evident once we differentiate (pTI). Successively inserting ( pH]) into 
(|57|), (|55[ ) and ( p5D we obtain the (Q/q)'^ contribution to the last term in square brackets in 
(|2|) and (0). 

It remains to evaluate the term containing F^^{Qti minjQ, q}) in (^5|). We may rewrite this 
term as 



- mQ{Ql-q')\F.{QT.q) 



Q{q-Q) 



+ F^{QT,Q)Q{q~Q)Y 



(59) 



We immediately identify the term with 1 in the square brackets as one of the terms in (P6l). We 
are left with two contributions from (^), each containing 0(g — Q), which combine to become 

as f^T d'^q 



UJ JQ'^ nq^ 



F^iQT,q)-F^{QT,Q) 



as rQTd'qdFUQT,q),_fq^ 



io Jq'^ vrg^ 91ogg^ 



log 



Q 



2 / ' 

T' 



(60) 



where the last equality follows on integrating by parts. We have now deduced (^6]), and (pT]) , 
from (ESl). 



Appendix B : Numerical technique used to solve CCFM equation 

We briefly describe the numerical technique that we used to solve the integral equation, 
(|i5|), for the gluon distribution F{x,Qt,Q)- The starting point is the (double) expansion of 
the function f{x,QT,Q) = Q'tF{x,QtiQ) in terms of the Tchebyshev polynomials with their 
arguments being linear functions of the variables log(Q|i/A^) and \og{Q'^ /h?). To be precise 
we map the regions Qq < Q\ < Qj^ and fi^ < < Qp into the interval (-1, 1) introducing 
respectively the variables 

We then expand the function / in the polynomial form 

N 

fix, Qt, Q) = E a(rT)C,(r)^.(x) (62) 
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where the functions fij{x) are the values of /(x, Qt, Q) at the "nodes" Q'ip^ and specified by 



Qti ( Qf\ ' Q] f Q 



F 



QfQo \Qo J Qf^- V ^^ 

with Tfc defined by 



(63) 



rfc = cos{(A;-i)7r/iV}. (64) 
The functions Ckir) in (|6^ ) are obtained from the Tchebyshev polynomials Tn{T) as follows 

1 ^ 

Cfc(^) = H VnTn{r)Tn{Tk) (65) 

-^^^ n=l 

where vi = \ and f „ = 1 for n > 1. In this way we can achieve a good approximation to the Qt 
and Q dependence of / in terms of a modest number N of Tchebyshev polynomials. Typically 
we take = 10. 

We now substitute the expansion (|6^ ) into the CCFM equation ( ^5|) and obtain the following 
set of Volterra-type integral equations for fij (x) 



Ux) = /°(x) + / - E (66) 
•^^ ^ k,i=i 

where the driving term 



f^,i^)=Ql^F%X,QT^,QJ), (67) 

and the kernel 

Jo J TTQ ^Ixi 

e{Q'S,-Ql)A^s{z,q,QT^,^^')Ck{T) Ceir). (68) 

Here we have 

Qri = \Qt + ^(1 -z)? = Qli + (1 - zYq"" + 2(1 - z)qQTi coscj) (69) 



and the non-Sudakov form factor Ajv^ of (^). The set of equations of (|66D is solved for the 
fij{x) by iteration and the function f{x,QT,Q) then determined from (|62D . 



Appendix C : Different forms of the BFKL equation 

Here we show the equivalence of the forms of the BFKL equation given in ( [5l| ) and ([531). 
We start from (^Tj) and write it in the form of a two-dimensional integration over cPQ'j,. It 
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is convenient to introduce a regulator 5^ and to combine the two virtual contributions which 
contain F{z,Qt)- Then (|5TD, with = 0, gives 

F{x,Qt) = F°'^{x,Qt) + 

. y _ f'dz f d'Q'r 1 Q'tF{z,Qt) 1 

(70) 

We introduce an integration over the Feynman parameter A to express the virtual contribution 
in (ITOD in the form 



I = Q 



d'Q'T 1 



\Q't - Qt? + ^^][Qt + \Q't - Qt? + <^'] 



TT 



If we now make the substitution ^ 

Qt ~^ Qt "I" _|_ 

then we can easily perform the d'^Q'j^ and dX integrations and obtain 

I = log f #U 0(5Vg^) 



52 7rg2 



+ 0{6'/Q'r). (72) 



We now substitute the virtual contribution / of ([7^) back into ( [TOl) and change the variable of 
integration over the real emission contribution to q = — Q^. In this way we obtain the 
form of the BFKL equation shown in (^). 



References 

[1] L.V. Gribov, E.M. Levin and M.G. Ryskin, Phys. Rep. 100 (1983) 1. 

[2] For recent reviews of small x physics see e.g.: A.D. Martin, Acta. Phys. Polon. B25 (1994) 
265; A.D. Martin, Proc. Int. Conf. on HE Physics, July 1994, Glasgow; J. Kwiecinski, Proc. 
of the Workshop QCD94, July 1994, Montpellier, France; E. Laenen and E.M. Levin, 
Fermilab-Pub-93-350-T; S. Catani, invited talk at Les Rencontres de Physique de la Valee 
d'Aosta, "Results and Perspectives in Particle Physics" , La Thuile, March 1994, Florence 
preprint DFF 207/6/94. 

[3] J. Feltesse, Proc. Int. Conf. on HE Physics, July 1994, eds. P. J. Bussey and I. G. Knowles, 
Glasgow. 



18 



[4] See for example, J. C. Collins, D. E. Soper and G. Sterman in "Perturbative QCD", ed. 
A. H. Mueller, (World Scientific, Singapore, 1989), p. 1, and refs. therein. 

[5] G. Altarelli and G. Parisi, Nucl. Phys. B126 (1977) 298; E. Reya, Phys. Rep. 69 (1981) 
195; G. Altarelli, Phys. Rep. 81 (1982) 1. 

[6] Yu. L. Dokshitzer, Sov. Phys. JETP 46 (1977) 641. 

[7] S. Catani, M. Ciafaloni and F. Hautmann, Phys. Lett. B242 (1990) 97; Nucl. Phys. B366 
(1991) 657; J. C. Collins and R. K. Ellis, Nucl. Phys. B360 (1991) 3; S. Catani and F. 
Hautmann, Nucl. Phys. B427 (1994) 475. 

[8] E. A. Kuraev, L.N. Lipatov and V. Fadin, Zh. Eksp. Teor. Fiz. 72 (1977) 373 [Sov. Phys. 
JETP 45 (1977) 199]; Ya. Ya. Balitskij and L.N. Lipatov, Yad. Fiz. 28 (1978) 1597 [Sov. 
J. Nucl. Phys. 28 (1978) 822]; L.N. Lipatov, in "Perturbative QCD", ed. A.H. Mueller, 
(World Scientific, Singapore, 1989), p.411; J.B. Bronzan and R.L. Sugar, Phys. Rev. D17 
(1978) 585. 

[9] T. Jaroszewicz, Acta. Phys. Polon. Bll (1980) 965. 

[10] A. J. Askew, J. Kwiecihski, A. D. Martin and P. J. Sutton, Phys. Rev. D49 (1994) 4402; 
A. J. Askew et al., Phys. Lett. B325 (1994) 212. 

[11] M. Ciafaloni, Nucl. Phys. 296 (1988) 49. 

[12] S. Catani, F. Fiorani and G. Marchesini, Phys. Lett. B234 (1990) 339; Nucl. Phys. B336 
(1990) 18. 

[13] G. Marchesini, Proceedings of the Workshop "QCD at 200 TeV" , Ericc, Italy, June 1990, 
edited by L. Cifarelli and Yu.L. Dokshitzer, Plenum Press, New York, 1992, p. 183. 

[14] G. Marchesini, Proceedings of the QCD94 Workshop, July 1994, MontpeUier, France. 

[15] G. Marchesini, Univ. of Milan preprint, IFUM-486-FT (November 1994). 

[16] B.R. Webber, Nucl. Phys. B (Proceedings Supplement) 18C (1990) 38. 

[17] G. Marchesini and B.R. Webber, Nucl. Phys. B386 (1992) 215; B.R. Webber, Proceedings 
of the Workshop "Physics at HERA", DESY, Hamburg, Germany, October 1992, edited 
by W. Buchmiiller and G. Ingelman. 



19 



Figure Captions 

Fig. 1 Gluon emission with probability given by (^. 

Fig. 2 Multigluon emission in deep inelastic scattering of a proton. 

Fig. 3 Kinematic variables for multigluon emission. The distribution JF„ describing n gluon 
emission is given in terms of the distribution Tn-i for ri — 1 gluon emission in (^). 

Fig. 4 The integrated gluon xg{x, Q"^) calculated using from the solution F(a;, Qt, Q) of the 
CCFM equation (^Sf) for different choices of the lower cut-off = /i^ at each of five values 
of Here we take Ql = 10^ GeVl Recall that our solutions are obtained from a "flat" 
gluon input. 

Fig. 5 The dependence of the gluon distribution obtained by solving the CCFM equation 
p5| ) for different choices of the lower cut-off = /i^ at each of two values of and four 
values of small x. The values of fx^ are as in Fig. 4. Here we take = 10^ GeV^. 

Fig. 6 The dependence of the gluon distribution F{x, Qt, Q), obtained by solving the CCFM 
equation (^) (continuous curves), compared with that from the DLL equation (dashed 
curves), for (a) Q'^ = 10 GeV^ and (b) Qt = 100 GeV^. In each case we show the 
dependence for x = 10~^, 10""^, 10~^ and 10~^. 

Fig. 7 The dependence of the gluon distribution F{x, Qt, Q), obtained by solving the CCFM 
equation ( ^5]) (continuous curves) and the DLL equation (dashed curves), for (a) Q^ = 10 
GeV^, (b) Q^ = 10^ GeV^ and (c) Q"^ = 10^ GeV^. In each case we show curves for 
X = 10-^ 10-^ 10-^ lO-^ and 10"^ 

Fig. 8 The integrated gluon distribution xg versus x, obtained from both the CCFM (continuous 
curves) and the DLL (dashed curves) integral equations, for Q^ = 4, 10, 10^, 10^ and 10'^ 
GeV^. Recall that our solutions are obtained from a "flat" gluon input. 

Fig. 9 The effective values of A, defined by xg = Ax~'^, obtained from the gluon distributions 
shown in Figs. 8 and 10. The CCFM values (continuous curves) are compared with those 
obtained from the BFKL (dot-dashed curves) and DLL approximations (dashed curves). 
In each case we show curves corresponding to the five different values of Q^ 

Fig. 10 The integrated gluon distribution xg versus x, obtained from the CCFM (continuous 
curves) and the BFKL (dashed curves) equations for = 4, 10, 10^, 10^ and 10^ GeV^. 
Recall that our solutions are obtained from a "fiat" gluon input. 
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